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Abstract 

We perform an unitary transformation for the symmetric phonon mode of the Holstein 
molecular crystal hamiltonian. We show how to compute the electronic spectral functions 
by exact numerical diagonalisation of an effective hamiltonian fully taking account of the 
symmetric phonon mode, usually discarded. 



This paper explains carefully how to simplify numerical investigations of the Holstein model, 
and other related model containing electron-phonon interactions, using exact diagonalization. 
Such a simplification was introduced many years ago by Ranninger and Thibblin for the two sites 
polaron problem]]]]. The symmetric phonon mode is picked out and taken account analytically 
as displaced oscillator, reducing the number of phonon modes by one (one is a large number 
for small clusters). Since then, this trick has been widely used. However, many authors^], 
neglected important contributions from the symmetric mode. 

We consider the Holstein hamiltonien, without Coulomb repulsion, given by 

J> j,$,<? j j,cr 
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The sum over 5 is on the z — 2d nearest neighbours (z — 1 for a two sites lattice) and 
e is some parameter. For a M sites lattice, the symmetric phonon mode is given by a s = 
(ai + a 2 + . . . + a M )/\/M and the corresponding hamiltonian is given by 

# s = w a;a s =- (a; + a a ), (2) 



where TV is the total number of electrons. The hamiltonian becomes H = H e + H s . H e is 
an effective hamiltonian with one phonon mode missing which has to be numerically solved. 
Indroducing new bosonic operators = a s — gN/y/M, we obtain the diagonal form 

H s = u b j N b N - g 2 N 2 u /M. (3) 

Let | n, N) be the eigenstates of H s with eigenvalues E^ N = mj — g 2 N 2 uj /M. This transfor- 
mation is an unitary transformation, bjq = a s with 

U N = e 9N(at-a s )/VM _ ^ 

The coherent states are given by | n, N) = Un \ n > where | n > are the eigenstates of aja s . 
The matrix elements are given by (n', N' \ n,N) — < n' \ Un-n 1 \ n >. In the particular case 
of a transition from the ground state with iV electrons to an excited state with N' electrons, 
we get 



We now show how this matrix element arise in the expression of the electronic correlation 
functions, say J\{t) = < c k)CT (t) c\. > where k is the momentum. Let | m, N > be the 
eigenstates of H e with eigenvalues E% N . At T = 0, we get 

.hit) =< 0, N | (0, N | c k , a (t) | 0, N > | 0, TV) 

= ^ < 0,7V | (0,7V | c Ka {t) | m, N + 1 > | n, N + 1) x 

m,n 

x < m, TV + 1 | (n, N + 1 | c[ ;(T | 0, TV > | 0, TV) (6) 
and the spectral function is given by 

Ji(u) = 2ttJ2 \(n,N + l\ 0,N) \ 2 \ < 0, N \ c Ka \ m,N + 1 >| 2 x 



x 8{u + E S ' N + E e ' N - E\ 
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This result is to be compared with the naive expression based only on H e which is 



J? aive (u,) = 2tt £ I < 0, N I c Ka | m, N + 1 >| 2 + ^ - E^ +1 ). 
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Figure 1 shows J\(k = 0,ou) with uj^/t = 0.2 and A = g 2 u>o/t = 1.2 for the two sites problem 
with ground state corresponding to iV = 0. Part (a) reproduces Alexandrov et als result (figure 
4. a of reference]^]) using equation (Ef), while part (b) shows the correct result using equation 



with E p = g 2 uj , an useful test formula. 

I thank J. Ranninger for interesting discussions. 
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Figure Captions 

Figure 1: Spectral function J\{k = 0,u) for the two sites Holstein Hamiltonian with 
uo/t = 0.2 and A = 1.2. (a) shows the naive result obtained without the symmetric contribution 
while (b) shows the correct result. 
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